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<N : ABSTRACT 

*^\ We present our latest results on near- to mid- infrared observation of SN2006jc 

O I at 200 days after the discovery using the Infrared Camera (IRC) on board 

^ \ AKARI. The near-infrared (2-5/im) spectrum of SN2006jc is obtained for the 

first time and is found to be well interpreted in terms of the thermal emission 
from amorphous carbon of 800±10K with the mass of 6.9 ± 0.5 x 1O~^M0 that 
was formed in the supernova ejecta. This dust mass newly formed in the ejecta 
Q I of SN 2006JC is in a range similar to those obtained for other several dust forming 

QQ I core collapse supernovae based on recent observations (i.e., 1O~^-1O~^M0). Mid- 

infrared photometric data with AKARI/IRC MIR-S/S7, S9W, and Sll bands 
have shown excess emission over the thermal emission by hot amorphous car- 
bon of 800K. This mid-infrared excess emission is likely to be accounted for by 
the emission from warm amorphous carbon dust of 320±10K with the mass of 
2.71q'5 ^ 1O~^M0 rather than by the band emission of astronomical silicate and/or 
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silica grains. This warm amorphous carbon dust is expected to have been formed 
in the mass loss wind associated with the Wolf-Rayet stellar activity before the 
SN explosion. Our result suggests that a significant amount of dust is condensed 
in the mass loss wind prior to the SN explosion. A possible contribution of 
emission bands by precursory SiO molecules in 7.5-9. 5/im is also suggested. 

Subject headings: dust, extinction - infrared: ISM - supernovae: general — 
supernovae: individual(SN2006jc) — stars: Wolf-Rayet 



INTRODUCTION 



Study on the dust formation in the ejecta of core-collapse supernovae (SNe) is an impor- 
tant topic to explore the origin of dust in the early universe. Besid es the theoretical studies 

that s uggest the dust condens ation in the ejecta of core-collapse SNe (JKozasa. Hasegawa &: Nomoto 
199ll : iTodini fc Ferrara]l200ll ). several pieces of the observational evidence for the dust for- 
mation in SN ejecta have so far been reported. For example, up to O.O2M0 of dust formation, 
which is close to the value of O.IMq n eeded for core-collapse sup ernovae to account for the 



dust content of high redshift galaxies (IMorgan fc Edmunds 



II supernova 2003gd has been reported by ISugerman et al. 



20031). in the ejecta of the type- 
J2006h . As for the case of the 



type-II supernova 1987A, however, even the highest estimate of the condensed mass in the 
SN ejecta re a ches only up to 7.5x10~'^Mq (JErcolano. Barlow fc Storeyl l2005l ). Moreover, 
Meikle et al.l (120071 ) revised the mass of newly condensed dust in the ejecta of SN 2003gd 
as 4xlO~^M0 based on recent observations with Spitzer. Spitzer Miltiband Imaging Pho- 
tometer (MIPS) observations of the Galactic core-collapse supernova remanant Gas A show 
that the dust mas s associated with it is much smaller (~ O.OO3M0) than previously thought 



Hines et a 



Rho et al. 



200J), while a recent Spitzer Inirared Spectrofraph (IRS) mapping observation 
20071 ) reports that the total mass formed in Gas A should be at least O.O2M0. 



There still remains a gap in the produced dust mass in core-c ollapse SN ejecta be tween those 
observational results and theoretical prediction of O.I-IM0 (JNozawa et al.ll2003l ). 



Supernova (SN) 2006jc is a pec uliar Type lb supernova and was discovered on 2006 Octo- 
ber 9.75 (UT) (JNakano et al.ll2006l ). It is believed that the progenitor star had experienced a 
luminous outburst s imilar to those of luminous blue var iables (LBVs) 2 years prior to the su- 
pernova event fiPastorello et al.l 120071 : iFoley et al.l 120071 ). Evidence of substantial interaction 
of the SN ejecta with a dense He-rich circumstellar medium (GS M) ejected during the LBV- 
like e ruption was reported based on Chandra X-ray observations (llmmler. Modjaz fc Brown 
2006h . 
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Recent observational studies on this supernova have shown that SN2006jc is an in- 
teresting target for the study of dust forming n iassive-star supernovae . The near-infrared 
re-brightening of SN2006JC was firstly reported by lArkharov et al.l (120061 ) from late November 
through early December 2 06 (~ 50 days after the discovery) based on J, if, and K band 
observations. ISmith et al.l (120071 ) report that both the appearance of a strong continuum 
emission at red/near-infrared wavelengths and the fading of redshifted sides of the narrow 
He I emission lines occurred simultaneously between 51 and 75 days after the brightness 
peak. These characteristics are interpreted as the evidence for dust formation in SN2006jc, 
although the timescale is much shorter than the general dust formation timescale, at least 
a few hundred days, typical of other dust forming SNe. They found that graphite grains 
with T ~ 1600K or slightly hotter silicate grains well fit the optical(red) spectrum of day 75, 
pointing out that the dust is mainly carbon and not silicate because of its high temperature. 

While most of the dust forming supernovae t hat have been observa tionally reported 
are type-II except for the case of type- lb SN 19901 (JElmhamdi et al.ll2004l ). the type-lb SN 
2006JC gives us a unique opportunity to investigate the dust formation not only in the SN 
ejecta but also in rich circumstellar materials that have come from the mass loss events 
prior to final core-collapse. Near- to mid-infrared observations of SN2006jc in the early 
phase of dust condensation are essential to examine the composition and properties of the 
newly formed dust in the SN ejecta and the pre-existing dust in the circumstellar medium 
separately. In this paper we present the new near-infrared spectrum and the mid-infra red 
photometric da ta of SN2006ic taken w ith the Infrared Camera(IRC) (jOnaka et al.ll2007l ) on 
board AKARI jMurakami et al.l 120071 ) on 29 April 2006 (200 days after the discovery). We 
derive the properties of the carriers of re-brightened near-infrared emission and the mid- 
infrared emission to understand the dust formation in SN2006JC. 
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2. Observation and Data Reduction 

Two pointed observations of SN2006jc were performed at 00:36:22(UT) and 02:15:47(UT) 
on 29 April 2007 as part of tlie director's time of AKARI. The former observation (ID:5124071) 
was performed with the AOT04 mode, in which the mid-infrared spectroscopic data were 
taken with two grisms, SGI (5.4-8.4yum) and SG2 (7.5-12. 9/im), and the mid-infrared imag- 
ing data were taken with the broad band fiher, S9W(6.7/im-11.6/im) of the MIR-S channel 
and the near-infrared spectroscopic data were taken with the grism, NG (2.5-5.0/im), and the 
near-infrared imaging data were taken with the broad band filter, N3(2.7-3.8/im) of the NIR 
channel. The latter observation (ID:512472) was performed with AOT02b mode, in which 
the mid-infrared imaging data were taken with two medium band filters, S7 (5.9-8.4/im) and 
Sll (8.5-11.3/im) of the MIR-S channel and simultaneously the near-infrared spectroscopic 
data was taken with the prism, NP (1.8-5.2/im), and the near-infrared imaging data were 
taken with the broad band filter, N3(2.7-3.8/im) of the NIR channel. The target was very 
faint and the spectrum was seriously affected by the blending with that of the host galaxy 
UGC4904, which hinders us to obtain useful spectral information with NG, SGI and SG2. 
In the following analysis we use the spectrum taken with NP as well as the imaging data 
with N3, S7, S9W, and Sll. The total exposure time was 206s, 206s, 224s, 56s, and 224s for 
NP, N3, S7, S9W and Sll, respectively. 

Each of imaging data reduction procedures of subtraction of the detector dark current, 
corrections for the effect of the hit by high-energy ionizing particles and the scattered light 



(jSakon et al.l 120071 ) . the fiat fielding, the distortion correction, and the shift and co-addition 
of the exposure frames for N3, S7, S9W, and Sll bands follows those in the IRC Imaging 
Pipeline Version 20070912 and the unit conversion factors from Analogue-to-Digital Con- 
verter Unit (ADU) per unit time to Jy are taken from AKARI IRC Data Users Manual 
ver.1.3. The pixel scale of the NIR channel is l."46 and that of the MIR-S channel is 
2." 34. The FWHM of the image size is 4."0, 5."1, 5." 5, and 4." 8 and the peak central pixel 
fiux is 5.9%, 12.0% , 11.5%, and 12.1% for NIR/N3, MIR-S/S7, S9W, and Sll, respectively 



(jOnaka et al.l 120071 ). Each of the spectroscopic data reduction procedures of subtraction of 
the detector dark current, correction for the high-energy ionizing particles effects, flat fleld- 
ing, and the shift and co-addition of the ex posure frames for NP data follows those in the 



IRC Spectroscopy Toolkit Version 20070913 flOhvama et al.l 120071 ) 



Fig. [Bi shows the RGB false color image of SN2006jc produced with AKARI/IRC 
NIR/N3(blue), MIR-S/S7(green), and Sll (red) bands, where the data taken with N3 and 
Sll bands are degraded into the gaussian beam with the FWHM of 5."1 so that they match 
with the image size of the S7 band data. SN2006jc is located close (~10") to the nucleus of 
the host galaxy UGC4904 and the diffuse background or foreground component of UGC4904 
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Fig. 1.— (a) RGB false color image of SN2006jc obtained with the AKARI/IRC NIR/N3 
(blue), MIR-S/S7 (green), and Sll (red) bands, (b) Slitless spectroscopic image with 
NIR/NP. The dispersion direction to longer wavelengths is shown with the white arrow. 



must be carefully removed to obtain the flux solely from the SN2006jc. The photometric 
decomposition technique employed is described in §2.1. The near infrared spectral image 
of SN2006JC is obtained by the slitless spectroscopy with AKARI/IRC NP (see Fig. lb) and 
it shows that the spectrum of SN2006jc is seriously blended with that of UGC4904 and a 
careful subtraction of the UGC4904 component is needed to derive the pure spectrum of 
SN2006JC. The spectroscopic decomposition technique employed is described also in §2.2. 



2.1. photometric decomposition technique 

To subtract the host galaxy component, we first derive the most appropriate point 
spread function (PSF) in the image and assume the galaxy component as a convolution of 
the source intensity distribution with PSF. This procedure is carried out in the vertical and 
horizontal direcion independently to check the reliability and consistensy. 

We obtain the intensity profile Pi{X) {i = 1, 2, ...n) in units of ADU pixel^^ along the 
i-th pixel row (see Fig. Wp), where n is defined by the number of the pixel rows that cover 
the SN2006JC (see Fig. [2^). We also define the reference pixel row which goes across only 
the host galaxy UGC4904 and obtain the intensity profile P{X) in units of ADU pixel"^ 
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Fig. 2. — The photometric decomposition techniques adopted for observation of SN2006jc 
with AKARI/mC NIR/N3 band. 
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along the reference pixel row (see Fig. [2)d). The position of the reference pixel row is selected 
so that P{X) has as sharp a profile as possible and, therefore, we use P{X) as the PSF. 
Then, we assume the intensity profile of the UGC4904 component p^GCAdOA^-^-^ ^^ Pi{X) as 
a convolution with the PSF written as 

oo 

pUCC.90.^X)= J2r^,ik)PiX-k), (1) 

fc=— oo 

where the distribution function rii{k) is assumed to have the gaussian form given by 

r]i{k) = hiexpi^-('^^j\. (2) 

The best fit parameters of hi, 6i, and o", {i = l,2,...,ra) are obtained so that the sum of 
\Pi{X) - pUGCi904(^x)\'^ for 50 < X < 69 and 79 < X < 90 becomes minimum, where the 
data of Pi{X) in 70 < X < 78 have contribution from the signal of SN2006jc and excluded. 
The intensity profile of the SN2006jc component P^ '''^(X) is, then, calculated as 

pSN2006jc^^^ = P,(X) - i^^f^t74904(^)_ (3) 

Finally, the total N3 band flux estimated from the analysis on the intensity profiles along 
the X direction, f^{N3), is derived as 



n „ 



(4) 



where a is the unit conversion factor from ADU to Jansky. 



In order to check the validity of the estimated flux density of SN2006jc using the intensity 
profiles along the X direction, the same decomposition technique is applied along the Y 
direction. We define m of target pixel columns so that they cover over the whole SN2006jc 
(see Fig. [2]i) and obtain the intensity profile Qj{Y) (j = 1, 2, ..., m) in units of ADU pixel"^ 
along the j-th pixel column (see Fig. [2]f). We assume the intensity profile of the UGC4904 
component QUGC4904(y^ j^ q,(j-^ ^g 



QUGC.90.^y^^ J2^^^k)QiY-k), (5) 

fc=— oo 

where Q{Y) is the PSF in the y direction represented by the intensity profile along the 
reference pixel column which is chosen in the same way as in the x direction (see Fig. [2^) 
and the distribution function rij{k) is assumed to have the gaussian profile. It is obtained so 



that the sum of \Qj{Y) -Qf^^^^'^^{Y)\'^ for 1 < F < 18 and 28 < F < 40 become minimum, 
where the data of Qj{Y) in 19 < y < 27 are contributed by the signal from SN2006jc. Then 
the total N3 band flux estimated from the analysis on the intensity profiles along the Y 
direction /J(A^3) is derived as 



fuim = a 



m „ 



(6) 



The obtained N3 band fluxes f^{N3) and f^{N3) are listed in Table. [H These two 
values are in good agreement within uncertainties. We adopted the weighted average of the 
two values taking account of the error as the final N3 band flux of SN2006jc. 

We obtain the flux densities of SN2006jc with MIR-S/S7, S9W, and Sll bands using 
the same decomposition technique as that adopted for NIR/N3 band, and the results are 
summarized in Table. [H The flux densities estimated from the intensity profiles along the X 
direction and Y direction are in good agreement within uncertainties for any of MIR-S/S7, 
S9W, and Sll bands. 

Table 1: Decomposed flux density of SN2006jc on 29 April 2007 (200 days after the discovery). 



band 


f^ 


/J 




(mJy) 


(mJy) 


NIR/N3 


0.386±0.046 


0.376±0.042 


MIR-S/S7 


0.586±0.063 


0.548±0.070 


MIR-S/S9W 


0.694±0.098 


0.735±0.105 


MIR-S/Sll 


0.495±0.061 


0.501±0.098 
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2.2. spectral decomposition technique 

The near infrared spectrum of SN2006jc is obtained by the slit-less spectroscopy with 
NIR/NP, in which the spectrum is dispersed in the X direction and a spectrum of a certain 
source is contaminated by the spectra of other sources aligning in the dispersion direction. 
Therefore, the spectrum of SN2006jc suffers severe blending with the light from UGC4904. 
We have developed a spectral decomposition technique for the spectral image of NIR/NP. 
We use the imaging data taken with NIR/N3 to investigate the structure of UGC4904 and 
to reproduce the spectral components of UGC4904 in the spectral image of NIR/NP, and 
have obtained a pure spectrum of SN2006jc. 

For the first step, the position in the Y direction of the NIR/N3 data is adjusted to that 
of the NIR/NP data with an accuracy of one-tenth of the pixel. 

We define the reference pixel row in the imaging data taken with NIR/N3 covering from 
Y=28 to Y=29 so that it goes across only UGC4904 and does not include SN2006jc (see 
Fig. [3^1). The intensity profile along the reference pixel row, /y[28:29](-^), is obtained by 
averaging the data from Y=28 to Y=29 (see Fig. [3]d). Then, we define the target pixel row 
covering from Y=20 to Y=25 so that it goes across the major part of SN2006jc (see Fig. ^) 
and the intensity profile along the target pixel row, /y[20:25](-^)5 is obtained by averaging the 
data from Y=20 to Y=25 (see Fig. [3t). 

We model the intensity profile of the UGC4904 component, lyS^^^^^^X), contained in 

^Y[20:25] (^) with 

00 

-^Y[20?25r(^) = 2J ^Y[20:25](^)-^Y[28:29](-^ -^), (7) 

fc=— 00 

where the free distribution function, ?7y[20:25] {k), is determined so that the sum of |/y[20:25] i^)^ 
Iy§o.25]^{X)\'^ for 50 < X < 69 and 79 < X < 100 becomes minimum, where the data of 
Pi{X) in 70 < X < 78 have contribution from the signal of SN2006jc and excluded (see 
Fig. ^). The obtained result of ^y[20:25](^) is shown in Fig. H^. 

Then we define the reference pixel row and the target pixel row in the spectral image of 
NIR/NP so that each has the same Y range as that defined in the image of NIR/N3. The 
spectral intensity profiles along the reference pixel row, 5'y[28:29](X) (see Fig. [3]F) and along 
the target pixel row, 5'y[20:25] {^) (see Fig. [3^) are obtained. 

Using the spectral intensity profile along the reference pixel row, 5'y[28:29](-^)5 and the 
obtained distribution function ?7y[20:25](^), the spectral intensity profile of the UGC4904 
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Fig. 3. — Tlie spectroscopic decomposition techniques adopted for observation of SN2006jc 
with AKARI/IRC NIR/NP band. 
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component, SyS^.2^^^{X), contained in 5'y[20:25] (-^) can be modeled with 



'S'y[20:25] i^) — 2-^ riY[20:2S\{k)SY[2S:2Oi\{X -k), 
k=—oo 



which is shown with the gray hne in Fig. [3^. The spectral intensity profile of SN2006jc, 
'S'Y[20:25f ''(^)' Contained in 5'y[20:25] (-^) is given by 



•^¥[20:25] ^^ > " '^Y[20:25]l-^ j " '^Y[20:25] \^ ) 



(9) 



and is shown with the thin black line in Fig. [3^. 
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Fig. 4. — The obtained distribution functions of (a) r/Y[20:25](^) for target pixel row and 
(b)?7Y[35:40](/i;) for tcst pixcl row. 

In order to evaluate the accuracy in the estimate of the spectral intensity profile of the 
UGC4904 component by this spectral decomposition technique, we define the test pixel row 
both in the image of NIR/N3 and NP covering from Y=35 to Y=40 so that it goes across 
the host galaxy UGC4904 (see Fig. [3^1 and e) and applied the same technique (see Figs. [3li 
and h). 

We model the intensity profile of the UGC4904 component, I^^^.l^^'^{X)^ contained in 

^Y[35:40] {X) with 



-^Y[35?40°^(^) = Z2 VY[35-A0]{k)lY[28:29]{X -k), 



(10) 



k=—oo 



where the free distribution function 77Y[35:4o](fc) is derived in the same way as ?7y[20:25](^) (see 
Fig. [311). The obtained result of ^y[35:40](^) is shown in Fig. |Dd. 
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The spectral intensity profile of the UGC4904 component, Syf^^.^Q^^i^X), is given by 

oo 
'5^[Sr(^) = E W[35.40](A:)5y[28:29](X - k) . (11) 



A:=— oo 



The obtained results of Sy^^.^^^'^i^X) and the residual profile given by 5'y[35:4o] (-'^) — 
Syi^^.^Q^'^iX) are shown with the gray line and the thin line in Fig. [3)i, respectively. The 
standard deviation of the residual profile from zero is 8.1(ADU/pixel), which confirms that 
the spectral intensity profile of the UGC4904 component is properly estimated by our spectral 
decomposition technique. The value of 8.1(ADU/pixel) is used to estimate the systematic 
error in the spectral decomposition procedure. 

We have made the wavelength calibration and the division by the system spectral re- 
sponse of NIR/NP for the obtained spectra l intensity pro f ile Sy ,^Q.^^,-''^(X) following the 



procedure in the IRC Spectroscopic toolkit (jOhyama et al.l 120071 ). The absolute flux cal- 



ibration for the NP spectrum of SN2006JC has been made by using the ratio of the total 
flux count of SN2006JC in /y[20:25](-^) to that obtained by the photometric decomposition 
technique in §2.1. 
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RESULTS 



3.1. The near- to mid- infrared characteristics of SN2006JC 



The results of the photometry of SN2006jc with AKARI/IRC N3, S7, S9W, and Sll 
bands are hsted in Table. [2j Near-infrared photometry o f SN2006ic with H and K bands 
was performed simultaneously on 28 and 29 April 2007 (JMinezakil lYoshiil ) using the mul- 
ticolor imaging photome ter (MIP) mounted on the MAGNUM 2m telescope at Haleakala 
Observatories in Hawaii (JKobayashi et al.lll998al Jbl). and the results are summarized also in 
Table. [1 



Table 2: Results of the photometry of SN2006jc on 29 April 2007 (200 days after the discov- 
ery). 



Instrument 


band 


Xhand (/Um) 


flux density (mJy) 


MAGNUM/MIP 


H 


1.6 


0.05±0.02 


MAGNUM/MIP 


K 


2.2 


0.14±0.02 


AKARI/IRC 


N3 


3.2 


0.38±0.04 


AKARI/IRC 


S7 


7.0 


0.57±0.07 


AKARI/IRC 


S9W 


9.0 


0.71±0.10 


AKARI/IRC 


Sll 


11.0 


0.50±0.10 



The obtained near-infrared spectrum of SN2006jc on day 200 together with the photo- 
metric results are plotted in Fig. [51 The near-infrared spectrum of SN2006jc is characterized 
by a continuum emission peaking around at ~4/im, which can be attributed to the thermal 
emission from dust grains. Possible atomic hydrogen recombination lines of Pf-^ at 2.87/im, 
Pf-e at 3.04/im, Pf-(5 at 3.30/im, Pf-7 at 3.74/im, Br-a at 4.07yum and Br-/5 at 2.63/im are 
recognized though they are significant only within 2-3o" (see Fig. [5]). We also note a small 
dent around at 4.6/im possibly attributed to the CO absorption band. 



3.2. Emission from the Newly-Formed Dust in the Near-Infrared 



Assuming that spherical dust grains of a certain kind X with a uniform particle radius 
ax and a total mass of Mx are located at a distance of R from the observer and that they 
emit optically thin thermal radiation of the temperature of Tx(K), the observed fiux density 
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Fig. 5. — The near-infrared spectrum of SN2006jc on day 200 after the discovery with 
AKARI/IRC NIR/NP. Photometric data taken with MAGNAM/MIP H, K bands (open 
square) and with AKARI/IRC NIR/N3, MIR-S/S7, S9W, Sll bands (sohd square) are 
plotted together. The l-a error level in the NP spectrum are also shown. 
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profile is given by 

-1 



f^W = Mx(^lnpWx^ 7rB^{X,Tx)QtW [^)\ (12) 

where Qx'^{X) is the absorption efficiency and px is the density of a dust particle of compo- 
sition X. The near-infrared spectral energy distribution (SED) of SN2006jc on day 200 of 
H and K band photometric data and the 2-5 /xm spectroscopic data is fitted with f^{X) for 
each case of the amorphous carbon (X = a. car.) and the astronomic al silicate (X = a.sil. 



as a composition of dust. We assume the distance of i? = 25.8Mpc (jPastorello et al.l 120071 ) 
for SN2006JC. 

As for the amourphous carbon case, the absorption efficiency, Q'^,car.W^ ^^^ ^ spherical 
amorphous c arbo n grain with a radius of aa.car. = 0.01/im is calculated from the optical 



constants of lEdol (119831 ). Assuming pa.car. = 2.26g cm~^, the ffi is carried out with the 
equilibrium temperature Ta.car. and the total mass Ma.car. being free parameters. The result 
of the ffi is shown in Fig. ^^ and the best ffi parameters of Ta.car. = 800 ± 10 K and 
Ma.car. = 6.9 ± 0.5 X IO'^Mq are obtained. 

The flux density at each AKARI/IRC imaging band is calculated for the model spectrum 
taking account of the color correction, which is given by 

fr'^'^iband) = '% , (13) 

where Rbandi^^) is the relative system spectral response of each AKARI/IKC band and 
^band = c/Xband correspouds to the reference frequency of e ach AKARI/IRC band defined 



with AAr3=3.2yum, As7=7.0/im, XsQw=Q-Opva., A5ii=11.0/im (lOnaka et al.l 120071 ). The model 
values that were converted into the AKARI/IRC calibration system of /"•™^'(A^3), /^•™^'(S'7), 
f^-'"''^iS9W), and f^-'"'''iSll) are plotted by crosses in Fig. EK. The observed data at S7, 
S9W, Sll bands have excess emission over /^•'='^'--(57), /^•™''-(59iy), and /;-™''-(5'll), respec- 
tively. The interpretation for the excess component is discussed in §4. 

As for the sihcate case, the absorption efficiency, Qa'^suXX), for a sp herical a s trono mical 



silicate grain with a radius of aa.car. = O.Ol/zm is taken from the values in iDraind (119851 ). and 
Pa.sii. = 3.3g cm~^ is assumed. The ffi was carried out in the same way as for the amorphous 
carbon grains. The result of the ffi is shown in Fig. Eb and the best fit parameters of 
Ta.sii. = 920 ± 10(K) and Ma.su. = 4.2 ± 0.3 x W'^Mq are obtained. The flux density at 
N3, S7, S9W, and Sll bands is calculated for the best ffi model spectrum taking account of 
the color correction and is shown in Fig. [6)d. The model values of /^•**'(S'7), /"•**'■ (5*91^), 
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and /"•"**'■ (S'll) largely exceed the observed flux density at S7, S9W, and Sll bands. The 
disagreement in the mid-infrared flux shows that astronomical silicate is not a likely carrier 
of the continuum in the 2-5/im region on day 200. 

We conclude that amorphous carbon is a likely carrier of the near-infrared continuum 
on day 200. The spectrum is well accounted for by the amorphous carbon dust with the 
temperature of 800±10 K and with the tota l mass of M„ r n.r. = 6 .9±0.5 x 10~^Mq. This result 



is consistent with the suggestion made by ISmith et al.l (120071 ) that the re-brightened near- 
infrared emission on 79days is likely to be carried by carbonaceous dust with an equilibrium 
temperature of 1600K, but not by silicate dust. The equilibrium temperature of 800K at on 
day 200 derived in our analysis is consistent with the scena rio that this compon ent is newly 
formed dust in the free-expanding ejecta of the supernova (JNozawa et al.l 120071 ). 
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Fig. 6. — Near-to mid-infrared spectral energy distribution (SED) of SN2006jc on day 200 
constructed by the H and K band photometric data(opera square) , the 2-5/xm spectroscopic 
data(j9/'us), and the mid- infrared photometric data.{solid square). The results of the best fit 
spectrum modeled with Eq. (TT2|) to reproduce the near-infrared data are shown for the case 
of (a) amorphous carbon (X = a. car.) and (b) astronomical silicate (X = a.sil.). The flux 
density at AKARI/IRC NIR/N3, MIR-S/S7, S9W, and Sll bands simulated for the model 
spectrum taking account of the color correction is shown with the thick crosses. 
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4. Identification of the Mid-Infrared Excess Component Carriers 

The mid- infrared excess component is seen over the model spectrum of the amorphous 
carbon of 800K at the S7, S9W and Sll bands. One possible candidate for the excess compo- 
nent is the dust with an emission band in the mid-infrared region carried by silicate and/or 
amorphous silica dust, which is discussed in §4.1. Another possible candidate for the excess 
component is the thermal emission from the dust with temperature lower than BOOK. The 
properties of the lower temperature dust are discussed in §4.2. In §5 we discuss which one 
of these two candidates is more likely. 

In the following analysis, we model the near- to mid-infrared spectrum as a combination 
of A^ kinds of dust components (X^; i = 1, ..., A^) including the amorphous carbon discussed 
in §3.2 and then the model spectrum is calculated as 

fr^'^X) = EMx. (I^PxA) nB^{X,TxJQtW {^)\ (14) 

where the fit parameters Mx^ and T^. are the mass and the temperature of the dust compo- 
nent Xi and Qx^{\), ax^, and pxi are the absorption efficiency, grain radius, and the density 
of dust Xi, respectively. 



4.1. Silicate and/or Amorphous Si02 Model 

One possible candidate for the excess emission in S7, S9W, and Sll bands is the dust 
band emission of the silicate and/or amorphous silica dust. In Fig. [7] are sh own t he ab - 



sorption efficiency profiles Q^'^(A) of a morphous car bon {X = a. car.) from lEdol (1l983l ) 



astro nomical silicat e {X = a.sil.) from iDraind (1l985l ). and amorphous Si02 {X = silica) 



from iPhilippI (119851 ) together with the system spectral response curve of the AKARI/IKC 
NIR/N3, MIR-S/S7, S9W, and Sll bands. While Q'^'^^ar.W shows a smooth fetureless profile 
without any band structures, Qa^suiX) and QluicaW show the band structure peaking at 
~9.5-10.0/im and ~8.5-9.0/im, respectively. Therefore, these band structures intrinsic to 
the silicate-related dust could contribute to the excess emission in S7, S9W, and Sll bands. 



Firstly, we assume a two component (A^ = 2) model of the amorphous carbon {Xi = 
a.car.) and the astronomical silicate {X2 = a.sil.) and the near- to mid-infrared data of 
SN2006JC on day 200 are fitted with Eq. (HI), where Ta.car. = 800K and Ma.car. = 6.9x IO^^Mq 
are fixed and Ta.su. and Ma.su. are taken as the free parameters. The grain radii of both com- 
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Fig. 7. — Absorption efficiency profile QxW ^^ (^) amorphous carbon (X = a. car.), 
(b) astronomical silicate {X = a.sil.), and (c) amorphous Si02 {X = silica), (d) The system 
spectral response curve of AKARI/IRC NIR/N3, MIR-S/S7, S9W, and Sll bands. 
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ponents are set as O.Olyum. The best fit model spectrum is shown in Fig. [8^ and we obtain 
Ta.sii. = 710 ± lOK and Ma.su. = 1.12 ± 0.02 x lO'^M^ 



/0. 



However, the flux densities in S7, S9W and Sll bands predicted for this two component 
(amorphous carbon + astronomical silicate) model taking account of the color correction 
cannot well reproduce the observed mid-infrared SED characterised by the enhanced S9W 
flux density (see Fig. ^). For one thing, the band structure of astronomical silicate peaking 
at ~9.5-10.0/im is almost fully included not only in S9W but also, sometimes more efficiently 
depending on the dust temperature, in Sll. For another, the predicted model spectrum is 
too low in 6 < A < 8.5/im, where the Sll band does not have its sensitivity but the S9W 
band does. Therefore, we need to have another dust component that has a band structure 
in A < 8.5/im in its absorption efficiency profile. Amorphous Si02 (silica) is one of the can- 
didates for such component (see Fig. [Tt). Moreover, silica should be able to be produced in 
SN ejecta (JNozawa et al.ll2003l . 120071 ). Its presence in the Cassiopeia A supernova remnant 
has been reported by a mid-infrared spectroscopic observation w ith Short Wavelength Spec- 
trometer (SWS) on board the Infrared Space O bservatory (ISO) (jDouvion. Lagage &: Pantin 
200ll ) and also by a recent S'prfzer observation (JRho et al.ll2007l ). 



We assume a two component (A^ = 2) model consisting of hot amorphous carbon [Xi = 
a. car.) and amorphous sihca (X2 = silica), where Ta,car. = 800K and Ma.car. = 6.9 x 1O~^M0 
are fixed and TsUica and MgUica are taken as the free parameters. The grain radii of both com- 
ponents are set to be O.Ol/zm and we assume the density of amorphous Si02 PsiUca = 2.62g 
cm~^. The best fit model spectrum is shown in Fig.[8)D and we obtain TsUica = 790±10K and 
Mgiiica = 6.5 ± 0.1 X 1O~^M0. This time, the predicted fiux densities are in good agreement 
with the mid-infrared SED characterized by the enhanced S9W fiux density. 
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Fig. 8. — The best fit model spectra with Eq. flT^ for the case of two component (A^ = 2) 
models made of (a) amorphous carbon {Xi = a. car.) and astronomical silicate {X2 = a.sil.) 
and (b) amorphous carbon (Xi = a. car) and amorphous Si02 {X2 = silica). Symbols are 
the same as in Fig. O 
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4.2. Two Temperature Amorphous Carbon Dust Model 

Another possible candidate for the excess component is the thermal emission from 
the featureless dust with temperature lower than 800K. Our near- to mid-infrared data 
of SN2006JC on day 200 are fitted with Eq. (IT^ assuming a two component (A = 2) model 
of hot amorphous carbon [Xi = h.a.car.) and warm amorphous carbon {X2 = w.a.car.). 
The temperature and the integrated mass of the former component are fixed to the values 
obtained in §3.2 and those of the latter component, T^.a.car. and M^.a.car., are taken as the 
free parameters. The same properties of the absorption efficiency profile, the grain radius 
and the density are assumed both for the warm and hot components. The fitting is made so 
that the fiux densities of the model spectrum at S7, S9W, and Sll bands that were converted 
into the AKARI/IRC calibration system should best reproduce the observed fiux density at 
S7, S9W, and Sll bands. The best fit model spectrum is shown in Fig. [9] and we obtain 
T^.a.car. = 320 ± lOK and M^,a.car. = 2.7+°;^ X lO'^M©. 

We note that the observed near- to mid-infrared SED has double peaked characteristics 
with one peak being located in the near-infrared and the other peak in the mid-infrared, and 
that the model fit with two-temperature amorphous carbon dust results in a considerable 
temperature gap between the two. This temperature gap can be explained if we assume 
the situation such that only the hot component with the temperature of 800K is the newly 
formed dust in the ejecta of SN2006jc and that the 320K warm component corresponds 
to thermal emission from the pre-existing circumstellar dust farther away from the newly 
formed dust with a dust-depleted region in between. The dust-depleted region may have 
been created by the dust evaporation due to the shock breakout in the initial phase of the 
SN explosion. 
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Fig. 9. — The results of the best fit spectrum with Eq. (fT^ are shown for the case of two 
component (A^ = 2) models made of two temperature amorphous carbons {Xi = h.a.car. 
and X2 = w.a.car.). Symbols are the same as in Fig. El 
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Discussion 



Observational evidence for the presense of silicate a nd/or slica grains formed in SN 



ejecta has been reported for the Gas A supernova remnant (JDouvion. Lagage &: Pantirul2001 



Rho et al.ll2007l . and references therei n) , and the presense o f silicate dust in the circumstellar 
medium is confirmed for SN 1987A (JBouchet et al.ll2006l ). The progenitor of type- lb SN 
2006JC is considered to be a Wolf-Rayet star, sufferin g the mass loss du ring its late evolution 
and forming a dense circumstellar medium rich in He (JFoley et al.ll2007 ) . The He layer of core 
collap se supernovae is suggested to contain more carbon than oxygen (JNomoto fc Hashimoto 
19881 ) and, esp ecially, the Wolf-Rayet star with strong mass l oss leads to thick C-rich CSM 
and envelope (JLimongi fc Chieffil l2006l : iTominaga et al.l 120071 ) . Therefore, the formation of 
silicate and /or silica grains in the C-rich circumstellar medium of SN2006jc is less likely (e.g., 
Nozawa et al.ll2003l ). 



Even if the silicate and/or silica dust is newly formed in the ejecta of SN2006jc, their 
temperature reaches only up to lOOK taking account of the optical absorption and th e 
collision with ejecta gas with lOOOK on day ~200 after the discovery (JNozawa et al.l 120071 ). 
AKARI/mC S7/S11 band colors calculated for astronomical silicate and silica with lOOK 
are /;-^*'(57)//^"-^^'-(511) = 5.2 x lO^^ and ff'^''{S7)/ff'''''{Sll) = 4.1 x lO'^, respectively, 
which are by far smaller than the S7/S11 band color of ~ 0.5 ±0.2 obtained for the observed 
mid- infrared excess component. Therefore, both the astronomical silicate and amorphous 
silica with lOOK cannot reproduce the flux density at S7, S9W, and Sll bands of the excess 
component and should not be the major carriers of the excess. Consequently, the mid- 
infrared excess component over the model spectrum of the amorphous carbon of 800K is 
likely to be IR light echo by amorphous carbon with 320±10 K in the circumstellar medium 
discussed in §4.2. However, the MIR photometric data with S7, S9W, and Sll bands are 
insensitive to the newly formed silicate of less than ~ IM© as long as the tempearute is 
~ 100 K. In this sense, the 6.9xlO~^M0 of amorphous carbon should be a lower limit of the 
newly formed dust mass in the ejecta of SN 2006jc. 

The formation scenario of silicate or silica dust in SN ejecta is illustrated as the fol- 
lowings; As the oxygen-rich ejecta expands and the gas temperature becomes lower than 
~ 3000K, silicon monoxide molecules are formed. Then, when the SiO molecular gas is 
cooled to 1200-1500 K, chemical reactions to condense silicate or silica dust become active 



( iNozawa et al.ll2003l ). Observational evidence for SiO m olecules being the p recursors of sili- 
cate or silica dust in supernova SN1987A is reported by iRoche et al.l (Il99ll ). in which they 
find the coincidence of the decrease in the mass of SiO molecules with the onset of the dust 



emission from the supernova. iKotak et al.l (120061 ) also report the presence of SiO molecules 
in a Type IIP supernova 2005af based on mid-infrared observations with the Spitzer Space 
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Telescope. 

We still find a small amount of excess in the S9W band over the best fit spectrum for 
the two temperature amorphous carbon dust model (see Fig. ^. SiO molecules carry the 
Av = 1 vibrational-rotational bands extending from 7.5-9.5/xm (JBeer. Lambert fc Sneden 



19741 ). which are more efficiently included in the S9W band than in the S7 and Sll bands. 
If the oxygen-rich part resides in the SN ejecta gas whose temperature is still higher than 
the condensation temperature of silicate and silica dust (~1200-1500 K), precursory SiO 
molecules can account for this excess in the S9W band. 

Finally, we conclude that amorphous carbon of 800±10K of 6.9±0.4x 10~^Mq is newly 
formed in the ejecta of SN2006jc on day 200. This dust mass is more than 3 orders of mag- 
nitude smaller than the amount ne eded for core collapse supe rnovae to contribute efficiently 
to the early-Universe dust budget (JMorgan fc Edmunds! l2003l ). Recent observational studies 
report the evidence of dust formation in the ejecta of several core collapse supernovae, but 
the produced dust mass in the SN e jecta is generally fou nd to be much smaller than theoret- 
ically predi cted values of 0.1-lM^ (iNozawa et al.ll2003l ). i.e. no more than 7.5xlO~^M(7) for 
SN 1 987A jErcolano. Barlow fc Storevlboosh. ~ IQ -^M^ for SN 1999em jElmhamdi et al. 
2QM), and AxIQ-^Mq for SN 2003gd jMeikle et al.ll2007h . The obtained amorphous carbon 



dust mass for SN2006jc on day 200 in our study is consistent with those relatively small 
values obtained for other several dust forming core collapse supernovae, although our data 
is almost insensitive to the newly formed silicate of ~100 K. It should also be noted that 
the dust mass is usually derived b y assuming that dust emission is optically thin, except for 



( jErcolano. Bar 



the dust mass ( Meikle et al. 



ow fc Storevl2005r) . and the derived mass is considered to be a lower hmit of 



20071 ). In addition to the dust formed in the ejecta, on the other 



hand, our mid-infrared photometric data suggest the presence of another amorphous carbon 
dust of 320 ± lOK of 2.7to5 ^ 10~^Mq as the circumstellar component. This component is 
expe cted to be formed in t he mass loss wind associated with the Wolf-Rayet stellar activity 
fe.g.. lWiUiams et al.lE992l : IWaters et al]ll997l : iMolster et aDll999l : IVoors et al.lboooh . It fol- 



lows that the dust condensation not only in the SN ejecta itself but also in the mass loss wind 
associated with the prior events to the SN explosion could make a significant contribution 
to the dust formation by a massive star in its whole evolutional history. 



6. Summary 

We present our latest results on near- to mid- infrared observations of supernova (SN) 
2006JC at 200 days after the discovery using the Infrared Camera (IRC) on board AKARI. 
The near- infrared (2-5/im) spectrum of SN 2006jc is obtained for the first time and is found 
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to be well interpreted in terms of the thermal emission from amorphous carbon of SOOdzlOK 
of 6.9 ± 0.5 X 10~^Mq that was formed in the SN ejecta. This newly formed dust mass 
in the SN ejecta is in a range similar to those obtained for other several dust forming core 
collapse SNe based on recent observations. Mid-infrared photometric data with AKARI/IKC 
MIR-S/S7, S9W, and Sll bands have shown excess emission over the thermal emission of 
amorphous carbon of 800K. This mid-infrared excess emission is likely to be accounted for 
by the emission from circumstellar amorphous carbon dust of 320±10K of 2.7]'^q'5 x 10~^Mq 
rather than by the band emission of astronomical silicate and/or silica grains. Since this 
circumstellar amorphous carbon dust of 320K is expected to be formed in the mass loss 
wind associated with the Wolf-Rayet stellar activity, our result suggests that a significant 
amount of dust is condensed in the mass loss wind associated with the prior events to the 
SN explosion. Finally, a possible contribution of emission bands by the precursory SiO 
molecules in 7.5-9.5/im is suggested to explain the enhanced S9W flux density. Mid-infrared 
spectroscopy of early-time SNe of similar type is highly important for further and detailed 
analysis on the properties of newly formed dust in the supernova ejecta distinguishing from 
the circumstellar or interstellar dust component. 

This work is based on observation of AKARI, a JAXA project with the participation of 
ESA. We thank all the members of the AKARI proiect, particularly those who have engaged 
in the observation planning and the satellite operation during the performance verification 
phase, for their continuous help and support. We would also express our gratitude to the 
AKARI data reduction team for their extensive work in developing data analysis pipelines. 
This work is supported in part by a Grant-in-Aid for Scientific Research on Priority Areas 
from the Ministry of Education, Culture, Sports, Science, and Technology of Japan and 
Grants-in-Aid for Scientific Research from the JSPS. 
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